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Abstract 

This work is devoted to the derivation of novel analytic results for special functions which are 
particularly useful in wireless communication theory. Capitalizing on recently reported series represen¬ 
tations for the Nuttall Q —function and the incomplete Toronto function, we derive closed-form upper 
bounds for the corresponding truncation error of these series as well as closed-form upper bounds that 
under certain cases become accurate approximations. The derived expressions are tight and their algebraic 
representation is rather convenient to handle analytically and numerically. Given that the Nuttall —Q and 
incomplete Toronto functions are not built-in in popular mathematical software packages, the proposed 
results are particularly useful in computing these functions when employed in applications relating to 
natural sciences and engineering, such as wireless communication over fading channels. 


I. Introduction 

It is widely known that special functions constitute invaluable mathematical tools in most fields 
of natural sciences and engineering. In the area of telecommunications, their utilization often 
allows the derivation of elegant closed-form representations for important performance measures 
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such as error probability, channel capacity and higher-order statistics (HOS). The computational 
realization of such expressions is typically straightforward since most special functions are built- 
in in popular mathematical software packages such as MAPLE, MATLAB and MATHEMATICA. 
Among others, the Marcum Q— function, Q m (a,b), the Nuttall Q— function, Q m , n (a,b ) and the 
incomplete Toronto function (ITF), T B (m,n,r) were proposed several decades ago and have 
been involved in analytic solutions of various scenarios in information and communication theory. 
More specifically, Marcum Q— function was firstly proposed by Marcum in [HI, 0 but it became 
widely known in digital communications, over fading or non-fading media, by the contributions in 
0-id. Its traditional and generalized form are denoted as Q\ (a, b) and Q m (a,b), respectively, 
and their properties and identities were presented in 0. Respective upper and lower bounds 
were reported in id- lUTl while an exponential integral representation was proposed in ifTTIl . In 
addition, useful closed-form representation for the Q m (a,b ) function for the special case that 
the order m is an positive odd multiple of 0.5, i.e. m ± \ G N, were derived in 0, [12]. 

Likewise, the Nuttall Q— function is a special function that emerges from the generalization of 
the Marcum Q— function 0. Its definition is given by a semi-infinite integral representation and 
it can be expressed in terms of Q m (a, b ) and the modified Bessel function of the first kind, I n (x), 
for the special case that the sum m + n is an odd integer lfl3l . Establishment of monotonicity 
criteria for Q m ,n{a, b) along with the derivation of tight lower and upper bounds and a closed- 
form expression for the case that m ± 0.5 G N and n ± 0.5 G N were reported in 0. In the 
same context, the incomplete Toronto function is a similar special function which was proposed 
by Hatley in lfl4ll . It is a generalization of the Toronto function, T(m. n. r), and includes the 
Q m (a,b) function as a special case. Its definition is also given in non-infinite integral form 
while alternative representations include two infinite series lfl5l . Its application has been used in 
studies relating to statistics, signal detection and estimation, radar systems and error probability 
analysis [fl6l -[H8l. 

Nevertheless, in spite of the evident importance of the Q m ,n(a,b ) and T B (m,n,r ) functions, 
they are not adequately tabulated while they are not included as standard built-in functions in 
popular software packages. Motivated by this, the Authors in 0, Il3l1l - ll34l derived explicit 
expressions for the Q m , n (a, b) and T B (m,n,r) functions. Specifically, a closed-form expression 
for Q m ,n(a,b) was derived in 0 for the specific case that m ± \ G N and m ± \ G N. 
A similar expression for the T B (m,nr) was subsequently derived in OT1 . If32]| . Furthermore, 
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simple series representations for both Qm,n(o>,b) and Tjfirri. n. r) functions were derived in 
If32l - lf34fl . These series are particularly useful since their algebraic form is relatively simple, 
which render the functions convenient to handle analytically. However, their form is infinite and 
thus the determination of adequate truncation needs to be taken into consideration for ensuring 
acceptable levels of accuracy for any given application. 

Motivated by this, this work is devoted to the derivation of novel upper bounds for the 
truncation error of the above functions. The derived representations are expressed in closed- 
form and have a relatively simple algebraic form. In addition, simple upper bounds are derived 
in closed-form, which are shown to be quite tight and in certain cases they become remarkably 
accurate closed-form approximations. 

II. Novel Results for the Nuttall Q-function 
A. Definition, Basic Properties and Existing Expressions 

The Nuttall Q-function is defined by the following semi-infinite integral representation (3j eq. 

( 86 )], 


Qm,n(a,b) = J x m e 2 I n (ax)dx, 

and constitutes a generalization of the Marcum Q-function, 


( 1 ) 


1 POO a 2 

Q m {a,b) = ——r x rn e~ 2L ^ ! ~I m _ l (ax)dx. (2) 

o m 1 J b 

The normalized Nuttall Q —function is given by Qm,n{a, b) = Q m ,n(a, b)/a n , which for n — 0 
yields Qig(a,b) = Qip(a,b) = Qi(a,b) = Q(a,b). For the special case that n — m — 1 it 
follows that Q m ,m-i(a,b) = Q m (a,b) and Q m , m _i(a, b) = a m ~ 1 Q m (a,b). Furthermore, when m 
and n are integers, the following recursion equation was reported in [13] eq. (3)], 


jjTn 1 I n (db') 

Qm,n(^5 a 2 +b 2 ^ ^ l)Qra—2,n(^) ^)* 


(3) 


e 2 


A finite series representation expressed in terms of the Qfia,b), I n {x) functions, the gamma 
function, r(x), and the generalized k th order Laguerre polynomial, L^‘\x), was proposed in 
[T3l eq. (8)]. However, this expression is restricted to the case that m + n is an odd positive 
integer. 







4 


As already mentioned in Section I, a simple expressions for the Q m , n (a, b ) function were 
reported recently by PHl . il34!l . This expression is given by, 


where 


and 


p a 21 e~^T{p + ty- 21 r ( m+n + 2 l+ 1 ,*p\ 

Qm,n b) — , n -m+ 2 i+l , : 77“ 

/=o +1 +1)2 2 r(p — i +1) 


/■ oo 

r(a) = / t a_1 exp(—t) dt 

Jo 


(4) 


(5) 


r»oo 

T(a, x) = / t a_1 exp(—t) dt 


( 6 ) 


denote the gamma and upper incomplete gamma functions, respectively. As, p —* oo, the 
terms T(p + l)p l ~ 21 /T(p — l + 1) vanish and © reduces to the following exact infinite series 
representation, 


°° a 21 e~^T ( m+n t 2l+ \ ^) 
Qm,n( a ’ b) = ^ 


TJ l\Y{n + l + l)2 "~ m 2 2 ' +1 ’ 


(7) 


For the special case that m, n £ Z + , the Ffa. .r;) function can be expressed in terms of the finite 
series representation in 11201 eq. (8.352.4)]. Hence, after necessary variable transformation © 
can be also expressed as, 




P L 


Aa 2l b 2k V{jp + l)p l ~ 2l T{^^- + l) 
l\k\T(n + l + l)T(p — l + l)2 /+fc ’ 


where, 


( 8 ) 


m + n — 1 
T = -~-1" 


(9) 


and 


^ = a”2 


m —n — 1 


(10) 


From ® we can also obtain straightforwardly, 
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oo L 


Aa?itfk V (m±z±L + l ) 

lvd^(^TTTl)2 [ + ^, 


771+71+1 


a„,„M) = EE 


(id 


Z=0 k =0 


It is noted here that the above representations have a convenient algebraic representation and are 
not restricted since they hold for all values of the involved parameters. However, an efficient way 
that will determine effectively after how many terms the above series can be truncated without 
loss of accuracy, is undoubtedly necessary. 

B. A Closed-Form Upper Bound for the Truncation Error 

The accuracy of © is proportional to the value of p and the series converges quickly. How¬ 
ever, deriving a convenient closed-form expression that is capable of determining the involved 
truncation error analytically is particularly advantageous. 

Lemma 1. For m,n,a £ R and b £ R + , the following closed-form upper bound for the 
truncation error holds, 




( 12 ) 


where, 



(13) 


with 



X 


t a 1 exp(—t) dt 


(14) 


denoting the lower incomplete gamma function. 
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Proof: The truncation error of © is expressed as, 

°° a 21 e~^T{p + l)p l ~ 21 T ( ™+n+2i+i ^ ^ 

e ' ~ , 5 i «r(n + i +1)2*=^ r(p - i + ij~ 

which can be equivalently expressed as 


(15) 


°° a 21 e~4r(p + Ijp 1 - 21 r f m + n + 2l +\ 

'< = E 

1=0 


tz /!r(n + ) + i)2 2 ^ ±i r(p-/ + i) 


Xi 


(16) 


p a 21 e-^Yf) + l)pl - 21 r ( m±n±2i±l ; 

S (!r(n + / + i)2 r, ~ m r +1 r(p-/ + i) 

As already mentioned, since the series in Z x tends to infinity, the terms T(jt + /)/; l2, /r(jt — / +1) 
vanish, yielding 


2-1 Qm,n(®) (j) ^ ^ 


°° a a e“^r ( m+n +2 ; +i ,f 


^ l\T{n + l+ l)2 n ^ ±1 ' 


(17) 


It is recalled that the normalized Nuttall O-function can be upper bounded by []9l eq. (19)], 
namely, 


Qm,n(0',b) < 2[" m lo.5,Mo.5 ( a > (18) 

where [n] 0 . 5 = \n — 0.5] + 0.5, with [.] denoting the integer ceiling function. To this effect, by 
substituting (fl8l) into (fTTl) and then into (fl6l) . one obtains, 


p a 21 e~^Y(p + l)p 1 - 21 T ( m+n+2l+ 1 ; ^ 
e t — ^Ho.5,Ho.s( a ^) _ “575 , , . r / , . 

i=o t!r(7i T- l + 1)2 2 T(p — l + 1) 


(19) 


Importantly, the upper bound in (IT8l) can be expressed in closed-form according to the expression 
in |[9l Corollary 1]. Therefore, after substitution in (fT9l) equation (fT2l) is deduced thus, completing 
the proollJ. ■ 


Remark 1. By following the same methodology as in Lemma 1, a similar upper bound can be 
deduced for the truncation error of the infinite series of Q m ,n( a jb) in ([XT]) . 

'By setting n = m — 1 in 0, 0, and t!2b . similar expressions are deduced for the Marcum Q-function, Q m {a, b). 
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C. A Tight Upper Bound and Approximation 

Besides the expressions for the special cases that m+n is an odd positive integer and m±0.5 £ 
N, n ± 0.5 £ N, no simple representations exist for the Nuttall Q-function. 

Proposition 1. For a, b,m,n £ M + and for the special cases that b —> 0 | a, m, n > § b, the 
following closed-form upper bound for the normalized Nuttall Q-function is valid, 


Qm,n (Oj b') f 




+ 1,4) 


n\ 2 2 e 2 


( 20 ) 


where \ F\ (a, b, x) denotes the Kummer confluent hypergeometric function 


Proof: According to Remark 2, equation © reduces to the exact infinite series as p —» oo. 
By also recalling that J 0 °° f(x)dx > / a °° f(x)dx when a £ K + , it follows that T(a, x) < r(x). 
As a result, the Q m ,n{a,b ) function can be straightforwardly upper bounded as follows: 


OO 


2=0 


2 / - 
a e 


p ^ m+n+2l+l ^ 


/!r(n + / + l)2 Ii ^ ±i 


“V* 

12 


( 21 ) 


By recalling that the Pochhammer symbol is defined as (a) n = T(a + n)/T(a) and expressing 
each gamma function as 


one obtains, 


r(x + l) = (xfTix) 


( 22 ) 


lo = 


T (^4f±±) 


e 2 


oo / m + n-\- I \ 

V 2 )l 


,21 


■ n — m + l 

n!2 2 


^ Z!(n+ 1),2* 

i=o v 


(23) 


The above infinite series can be expressed in terms of the Rummer’s confluent hypergeometric 
function 


iFi(a, b, x) = ^2 


( a)i x l 

to ^ ll 


Hence, by performing the necessary change of variables and substituting 
equation (l20l) is deduced and thus the proof is completed. 


(24) 
into (ED, 
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Fig. 1. Qm,n(a,b) in 0. 


Remark 2. When a , m, n > | b, the upper bound in (l20l) becomes an accurate closed-form 
approximation for Q,n : n((C b), namely, 






( m+yt+l 
2 


,n+ 1 , t ) 


I n — m-|-l q2 

n! 2 2 e 2 


(25) 


The Nuttall Q—function is neither tabulated nor built-in in popular mathematical software 
packages like MAPLE, MATLAB and MATHEMATICA. Hence, the derived expressions are 
rather useful both analytically and numericallj^. 

The behaviour of © is illustrated in Fig. 1 along with results obtained from numerical 
integrations. The series was truncated after 20 terms and one can notice the excellent agreement 
between analytical and numerical results. This is also verified through the level of the involved 


-li 


absolute relative error e r =| Q m , n (a, b) — Q m , n (a , b ) | / Q m ,n(a, b) which is less than e r < 10 
Similarly, the behaviour of (l20l) is depicted in Fig. 2. Clearly, it upper bounds the Q m , rt (a, b) 


2 Similar expressions for the Qm,n(a, b) function can be straightforwardly deduced with the aid of the identity Q ra , n (a, b) 
Qm,n{a,b) / a" i.e. by multiplying equations 0, 0, d!2t and ( 121b with a n . 















tightly while it becomes an accurate approximation for higher values of a. 



Fig. 2. Q m ,n(a,b) in ((20). 

III. Novel Analytic Results for the Incomplete Toronto Function 
A. Definition and Basic Properties 

The ITF has been also useful in telecommunications and is defined as, 

r B 

T B (m,n,r) = 2r n ~ m+1 e~ r2 / t m - n e~ t2 I n {2rt)dt. (26) 

J o 

When B = oo, the Tffirn, n. r) function reduces to the T(rn, n, r) function, while for the specific 
case n — (m — l)/2 it is expressed in terms of the Marcum Q-function, namely, 

T b ^ m , m ^ = 1 - Qm±i (rV 2, Bs/^j . (27) 

Alternative representations to the T B (a,b,r ) function include two series which are infinite and 
no study has been reported on their convergence and truncation [fT5ll . In addition, the following 
polynomial representation was proposed in ll32l . 
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p 

T B (m, n, r ) ~ 

fc =0 


r (p + k)r ^ n+k ^- m+1 7 p±i + k, B 2 ) 
/c!p 2fc -ir(p — k + l)T(n + k + l)e r2 


which as p —> oo it reduces to, 


T B (m, n, r ) = ^ 




r , 2( n+fe)-m+l 7 (m±l + ^ £2) 

k\T(n + k + l)e r2 


(28) 


(29) 


A Closed-Form Upper Bound for the Truncation Error 

A tight upper bound for the truncation error of (l29l) can be derived in closed-form. 
Lemma 2. For m,n,r el, Be M + and m> n the following closed-form inequality holds, 


„ [nJ “ a r <2t+,) ([njo.5 + fe— 5 )! (i — fc)! 

£ ‘“ fcj h *!« (L«Jo. 5 - fc - i)!(£ - fc - 0! 

7 [iU, (B + r) 2 ] 7 [!± 1 , (B - r) 2 ] 

^_l)M-^2 2fe+ i + (—l) fc 2 2fc+1 

^ T(p + fc)r 2 ( n+fc )- m+1 7 (^±1 + k, B 2 ) 

“ /c!p 2 fc- 1 r(p — k + l)r(n + k + l)e r2 

Proof: Since the corresponding truncation error is expressed as 

OO OO p 

u = ^2f(x) = J2f(x) -^2f(x) ( 31 ) 

p+1 1=0 l 

and given that (l29l) reduces to an exact infinite series as p —y 00 , it follows that, 




e t = 


E 

k =0 

7 _ 


j.2(n ■ k) m ■ 1 „ /m+1 _J_ k q2\ 

k\T(n + k + l)e r2 
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E 


k =0 


r(p + k)r 2 ( n+ V~ m+1 7 ( 2 ^i + k, B 2 ) 
k\p 2k ~ x Y {p — k + l)T(n + k + l)e r2 


Given that X 3 = T B (rn, n, r), the e t can be upper bounded as follows: 


(32) 


p 

e t < T B {\m\, H0.5,0 - 

k=0 


r(p + A;)r 2 ( n+fc )- m+1 7 + jfe, £ 2 ) 

^!p 2 fc-ip('p _ _|_ l)T(n + k + l)e r2 


(33) 
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The T b (|~m], |_//J 0 5 . r) function can be expressed in closed-form according to the closed-form 
expression in lf32ll . Hence, by substituting in (l33l) equation (l30l) is obtained, which completes 
the proof. ■ 

Remark 3. By omitting the terms Y(p + k)jY lk /Y(p — k + 1) in (l30l) . a closed-form upper 
bound can be deduced for the truncation error of the exact infinite series in Remark 2. 


C. A Tight Closed-form Upper Bound and Approximation 

Capitalizing on the algebraic form of the T B (m,n,r) function, a simple closed-form upper 
bound is derived which in certain cases becomes an accurate approximation. 

Proposition 2. For m, n, r G M, B e R + and m, n, r < -7, the following inequality holds, 


T B fn, n, r ) < 


r (^) iFi(^i,n+l,r 2 ) 


(34) 


r m-2n- lT( n + l)e r 

Proof: It is recalled that the 7 (a,x) function can be upper bounded by the T(a) function 


since 


T(a) = / t a 1 exp(— t)dt = / t a 1 exp(—t)dt = y(a,x = oo) 

Jo Jo 

Therefore, the T B (m,n,r) function can be upper bounded as, 


(35) 


~ r 2(n+fc)-m+l r fin±l +k \ 

T B (m,n,r) < > -— - , v 2 . — 5 —~ 

K ’ ’ ; ^ k\Y(n + k + l)e r 


(36) 


fc =0 


which with the aid of the identity 


can be expressed as, 


T(a,n) = (a) n T(a) 


(37) 


T B (m, n, r ) < 


^,2n—tn+lp ^ m +1 j 

T(n + 1 ) e r2 


OO 


E 


YF) t r* 


(■n + 1 )k k\ ' 


(38) 


Importantly, the above series can be expressed in closed-form in terms of the Kummer’s confluent 
hypergeometric function. Thus, equation (l34l) is deduced and the proof is completed. ■ 
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Fig. 3. TB{m,n,r) in (23) 


Remark 4. The proof can be also completed by assuming B —> oo in (l26l) and utilizing / l20l 
ec/. (8.406.3)] and H 201 eq. (6.631.1)]. The incomplete Toronto function reduces then to the 
Toronto function, as the two functions are related by the identity T B=00 (7n,n,r ) = T(m,7i,r). 
Furthermore (f34l) becomes an accurate approximation when m, n, r < 2 B, namely, 


T B (m, 77 ,, r ) ~ 


r 


(12±1) lFl (ni±l jn + iy) 

r m-2n-iY(n + l)e r2 


(39) 


Equation (l29l) is depicted in Fig.3 along with results from corresponding numerical integra¬ 
tions.The agreement between analytical and numerical results is excellent and the relative error 
for (|29|) is e r < 10 -4 for truncation after 20 terms. In the same context (28) is depicted in 
Fig.4 along with numerical results for three different scenarios. The involved relative error is 
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proportional to the value of r and is e r < 10 6 when r < 1. 



Fig. 4. Ts(m, n, r) in (28) 


To the best of the Authors’ knowledge, the offered results have not been reported in the 
open technical literature and can be particularly useful in emerging wireless technologies such 
as cognitive radio, cooperative communications, MIMO systems, digital communications over 
fading channels, ultasound and free-space-optical communications, among others lf3Tl - lt54i . and 
the references therein. 


IV. Conclusion 

Novel upper bounds were derived for the truncation error of representations for the Nuttall 
O-function and the incomplete Toronto function. These expressions are given in closed-form 
and are have a tractable algebraic form. Since the Nuttall —Q and incomplete toronto functions 
are not included as built-in functions in popular mathematical software packages, the offered 
results are useful in computing these functions efficiently. As a result, the accurate computation 
of critical performance measures in digital communications that involve these functions becomes 
more feasible. 
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